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Nonintegrability and fluctuation by symmetry violation 
in perturbed harmonic oscillator systems 
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General Arts, Hiroshima National College of Maritime Technology, J^212-l, 
Higashino, Oosakikamishima-cho, Toyota-gun, Hiroshima 725-0231, Japan 

The violation of symmetry between the time series of elongation and contraction rates of 
phase-space point spacings is studied to examine the chaos in perturbed harmonic oscillator 
systems. A transition from integrability to nonintegrability is quantitatively evaluated by 
introducing a correlation coefficient for the degree of symmetry violation. A feature of 
the correlation coefficient is discussed, suggesting a possible advantage over the Lyapunov 
exponent. The relationship between a fluctuation property based on symmetry violation and 
nonintegrability is also discussed. 

§1. Introduction 

We demonstrate the characterization of chaos by exploiting a correlation coef- 
ficient for a rate of elongations and contractions of consecutive phase-space point 
spacings. In the case of an integrable orbit on a Poincare section plane/)' it is 
found that a degree of symmetry occurs between the elongation and contraction 
rates of the consecutive phase-space point spacings. When this symmetry occurs, 
the absolute value of the correlation coefficient becomes 1.0. On the other hand, 
the symmetry of the rates of elongation and contraction is violated when integration 
becomes impossible, and the absolute values of correlation coefficients are close to 
zero. This analysis provides us with a new method of considering how the nonin- 
tegrabilities depend on the degree of symmetry violation between the elongations 
and contractions of consecutive phase-space point spacings in perturbed harmonic 
oscillator systems. We compare the information obtained from the absolute value 
of the correlation coefficient with that obtained from the Lyapunov exponent^)' 
and show the advantage of using the correlation coefficient. Moreover, we show that 
a fluctuation property derived from the consecutive phase-space point spacings fol- 
lows the shape of a non-Gaussian distribution function as the degree of symmetry 
violation increases. 



§2. Consecutive phase-space point spacings 

To begin with, let us define In- The spacing of consecutive phase-space points on 
the Poincare section plane is defined as /„ = \/ {qn+i — InY + {Pn+i ~ PnY- Here, 
n enumerates a set of points on the Poincare section plane. The scaled spacings of 
is defined as Ln = In ' {ln)~^, where the bar denotes a "spacing average" over 
all spacings. The spacings were scaled to make Ln = 1.0. The following modified 
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nondimensional SU{3) Hamiltonian^)'^) was adopted in analyzing L„: 



where r = 2Af — q^^^'^ — p^-*-)^ — 5(2)^ _ ^(2)^^ g^j^^ y(i) g^j^^j y(2) ^^^gj^Q^g ^j^g strength 
of perturbation interaction in each degree of freedom. The system with the Hamil- 
tonian in Eq. (2.1) consists of three single-particle levels with energies 
and £^(^), where each level has A/'-fold degeneracy. All subsequent analysis assumes 
that = 1/(2) ^ y_ 

our numerical calculation, we used A/" = 30, and the 

single-particle energies were 

£(0) = 0, £:(^) = 1, and £:(2) = 2. The position at time t 
is written as qn \ and the position at time t + e is written as Qnln where e represents 
the time step of our simulation. 

The succession of mappings with the Hamiltonian in Eq. (2.1) is approximated 

as 

(2) _ (2) , g-H [qh' ,Pii';qn' ,Ph' ) 

J2) _ „(2) _ — Vin ^■i" / (r,n^ 

Pn+1 — Pn (2) ^ • ^/ 

Hereafter the discussion focuses on g(^) and p^^^ , which are written as q and p for 
simplicity. 
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Fig. 1. Poincare section plane for the SU{2,) Hamiltonian with (a) V = -0.009, (b) V = -0.045, 
and (c) V = -0.110. 
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If the interaction V is sufficiently small, the Poincare section map turns out to be 
closed curves, as shown in Fig. 1(a). Many orbits form invariant curves in this map, 
like Fig. 1(a), as guaranteed by the KAM theorem. ^^'^^ In this weak-interaction case, 
the time dependence of L„ turns out to be periodic, as shown in Fig. 2(a) and Fig. 
3(a). The alternating behavior of L„ for this case is shown in more detail in Fig. 2(d). 
The value of L„ is distributed in the range 0.943 to 1.03, giving a distribution width 
of 0.087. As is clearly seen from Fig. 3(a), the Fourier components'^ expressing the 
mixing effects from other periodic orbits are small in the case of a small interaction 
V = —0.009. Comparing Fig. 3(a) with Fig. 3(b), one can see that the peak of 
the Fourier coefficients becomes continuous when the interaction strength increases 
from V = —0.009 to V = —0.045. As the interaction strength increases further, the 
chaotic region begins to spread out from the hyperbolic fixed points. This is shown 
in Fig. 1(c). In this chaotic case, the value of L„ is distributed in a wider range from 
0.0 to 2.0, giving a distribution width of 2.0, and the periodicity of L„ disappears, 
as shown in Fig. 2(c) and Fig. 3(c). A set of peaks extending from n = to n = 5018 
and beyond with no clear maxima in the Fourier coefficients of consecutive point 
spacings reflects the nonintegrability of the orbit. 



Fig. 2. Consecutive phase-space point spacings computed using (a) V — —0.009, (b) V — —0.045, 
and (c) V = —0.110. The data set of Fig. l(a)-(c) was used, (d) Short-time development of L„ 
for V = -0.009. 
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Fig. 3. Root-mean-square of Fourier cosine coefBcicnts and Fourier sine coefficients for the consec- 
utive phase-space point spacings depicted in Fig. 2. (a) V = —0.009, (b) V = —0.045, and (c) 
V = -0.110. 

§3. Symmetry violation for consecutive phase-space point spacings and 

the relationship to chaos 

Wc will now examine a quantitative estimate of chaos that uses the absolute 
value of the correlation coefficient for the rates of elongation and contraction of L„. 
In the case where ^^^±1 > 1 the rate of elongation is defined as Rm = and 

in the case where ^^P^ < 1, the rate of contraction is defined as Vm = -^^r^- For 
example, the values g = 1.1, g = 1.2, g = 0.9, g = 0.8, g = 1.3, and g = 0.7 
give rise to elongation rates Ri = 1.1, R2 = 1.2, and i?3 = 1.3 and contraction rates 
ri = 0.9, r2 = 0.8, and r-^ = 0.7. The average rates of elongation and contraction for 
a number M of phase-space point spacings are given by 

^ M 1 ^ 

^=mE^-'^"=mE^- (3-1) 

m=l m=l 

One may introduce a correlation coefficient^^ for R^ and that takes the form 

M 

Y,{Rm-R){rm-f) 

^=^= ■ (3-2) 

M 
m=l 

Now, let us begin our analysis by considering the regular-to-chaos transition visible 
in Figs. l(a)-(c) and Figs. 2(a)-(c). Since the consecutive phase-space point spacings 
are periodic in the case of a weak interaction, the time dependences of Rm and are 
symmetric with respect to ^^jj^ = 1, as shown in Fig. 4(a). In this case, the absolute 
value of the correlation coefficient is close to 1.0. In Fig. 4, note that the time series 
of Rm and are independent of each other. As the interaction strength increases 
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only slightly, the orbit appears as shown in Fig. 1(b), and Rm and now behave 
as shown in Fig. 4(b). This orbit is, however, not integrable and it can be seen that 
the symmetry between Rm and is approximately maintained as shown in Figs. 
4(b) and 4(d). When the interaction strength increases as shown in Fig. 1(c), one 
may not expect symmetric behavior between R^n and r^, and this is demonstrated 
in Fig. 4(c). In this case, as seen from Fig. 5, the absolute values of the correlation 
coefficients are close to zero. 



Fig. 4. Elongation rate and contraction rate of phase-space point spacings. (a) V — —0.009, (b) 
V = -0.045, and (c) V = -0.110. (d) Short-time development of for V = -0.045. 
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Fig. 5. Graph of both the Lyapunov exponent and the absolute value of the correlation coefficient 
in the transition process from integrability to nonintegrability. Details of the comparison are 
described in section 4. 



§4. Features of the correlation coefficient analysis 

In this section, the process is analyzed using the absolute value of the correlation 
coefficient and the Lyapunov exponent, for the case in which the periodic orbit 
transitions to the nonintcgrablc orbit. 

In the case with weak interaction, the absolute value of the correlation coefficient 
quickly reaches 1.0 as shown in Fig. 6(a). When the interaction strength increases, 
the absolute values of the correlation coefficients are close to 0.0, as shown in the 
plot for V = -0.110 in Fig. 6(a). 



Fig. 6. (a) Time dependence of the absolute value of the correlation coefficient for Figs. l(a)-(c). 
(b) Time dependence of the Lyapunov exponent for Figs. l(a)-(c). The number of data points 
in (a) is roughly half that in (b) duo to the division into elongation and contraction rates used 

to calculate the correlation coefficient. 

When the orbit is completely integrable, the absolute value of the correlation 
coefficient has exactly 1.0 as an upper bound. On the other hand, if the orbit is 
nonintegrable, it has ICi^^r^l = 0.0 as a lower bound. The Lyapunov exponent 
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takes negative values when the system is integrable and positive values when the 
system is non-integrable. However, the upper and lower limits of the value of the 
Lyapunov exponent are not clear, as depicted in Fig. 6(b). 

§5. Asymmetric behavior in consecutive phase-space point spacings and 
nonintegrability in the Henon-Heiles system 

To date, the analysis of IC/?^,.^! has been performed using the specific SU{3) 
Hamiltonian with fixed energy £'=40. To confirm the usefulness of our analysis using 
\^Rmrm\j let us apply our method to the Henon-Heiles system^) and also study what 
happens when the energy of the system is varied. 

The Hamiltonian is given by 

(5-1) 

If the energy E is sufficiently small, the Poincare section map turns out to be 
integrable, as shown in Fig. 7(a). In this case, it is found that a degree of symmetry 
occurs between the Rm and as shown in Fig. 8(a). Since this symmetry occurs, 
the absolute value of the correlation coefficient becomes 0.5091691 shown in Fig. 
9(a). In Fig. 7(b), the successive intersections of the trajectory with the Poincare 
section plane lie on twenteen closed curves with ten hyperbolic fixed points located 
at the joins of adjacent closed curves. Nonintegrability appears in the Hcnon-Hcilcs 
system as energy increases. An example of this at = 0.1669 is shown in Fig. 
7(c). From Figs. 8(c) and 9(a), it can be seen that the symmetry between and 
rm is violated with the result that the absolute value of the correlation coefficient 
approaches zero. 



Fig. 7. Surfaces of sections for the Henon-Heiles Hamiltonian with (a) E = 0.005, (b) E = 
0.11796899, and (c) E = 0.1669 from numerical integration. 

Next, we describe an advantage that the correlation coefficient has when com- 
pared with the Lyapunov exponent. The unstable orbit of Fig. 7(b) is divided into 
two orbits, as depicted in Figs. 10(a)-(b). The absolute value of the correlation 
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Fig. 8. Elongation rate and contraction rate of phase-space point spacings. (a) E = 0.005, (b) 
E = 0.11796899, and (c) E = 0.1669. (d) Short-time development of for E = 0.11796899. 

coefficient for the orbit in Fig. 7(b) is represented by the short dashed line corre- 
sponding to E = 0.11796899 in Fig. 9(a). The value of \CR^rm\ markedly varies 
in the region around m 2450 as seen in Fig. 9(a). The corresponding values for 
the Lyapunov exponent arc represented by the short dashed line in Fig. 9(b). As 
can be seen in those figures, there is a clear change in the absolute value of the 
correlation coefficient between the sections corresponding to the orbit of Fig. 10(a) 
and the orbit of Fig. 10(b). The Lyapunov exponent calculation shows the value 
does not vary markedly with increasing time, as illustrated in Fig. 9(b). In this 
case, the differing orbits of the two time sections shown in Figs. 10(a)-(b) cannot be 
distinguished clearly. This suggest a possible advantage of the correlation coefficient 
over the Lyapunov exponent. The numerical results discussed here show that it is 
possible to use |Cij^7-ml sensitively detect even slight changes of orbit, something 
that the Lyapunov exponent is unable to do. 
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Fig. 9. (a) Time dependence of the absolute value of the correlation coefficient corresponding to 
Figs. 8(a)-(c). (b) Calculation of the Lyapunov exponent as a function of n. 



Fig. 10. The unstable orbit as shown in Fig. 7(b) is composed of orbit (a), which exists in the 
interval n = to n = 4900, and orbit (b), which exists in the interval n = 4901 to n = 9990. 
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§6. Symmetry violation and probability density distributions 

The main topic of this section is the determination of a relationship between a 
fluctuating property measuring symmetry violation and nonintegrability. To begin 
with, we define a new variable Om = Rm x fm- When the perturbation interaction 
is sufficiently small, as in Fig. 4(a), this quantity is stationary with a value of 
Om ~ 1.0 as shown in Fig. 11(a). This corresponds to the case where there is 
symmetry between and r„j and IC/j^r-ml ~ l-O- This suggests the following 
equation for the time evolution m when the system is completely integrable. 



Fig. 11. Computation of the fluctuation property Om = Rm x rm for perturbation interaction 
values (a) V = -0.009, (b) V = -0.045 and (c) V = -0.110. 

In the case of Figs. 4(b)-(c), the balance between the elongation rate and con- 
traction rate of the phase-space point spacings is broken, and this is confirmed in 
Figs. ll(b)-(c) as a fluctuation where ^ 0. The statistical properties of this 
fluctuation are now discussed. To handle the logarithmic behavior between the elon- 
gation rate and contraction rates, and hence of Om, let us introduce the variable 
Qm, where 





(6-2) 
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Fig. 12. Probability density distribution of Qm- The width of the histogram dQm is 0.05. (a) 
V = -0.009, (b) V = -0.045, and (c) V = -0.110. The dotted line is the function with the 
paremeter a = 14.5 in Eq. (6.3). 



For example, the values Oi = 2.0, O2 = 5.0, O3 = 10.0, O4 = 0.5, O5 = 0.2, and 
Oq = 0.1 correspond to Qi = 1.0, Q2 = 4.0, Q3 = 9.0, Q4 = -1.0, Q5 = -4.0, 
and Qq = —9.0. This variable behaves in a manner similar to Om for symmetry 
violations. 

Wc now investigate the distribution of Q„i with increasing nonintegrability. Fig- 
ures. 12(a)-(c) show the probability density P{Qm) of finding for the components a 
value between and Qm + dQm for each interaction value. The probability density 
of Qm behaves as a delta function when the system is integrable, as shown in Fig. 
12(a). On the other hand, it becomes a non-Gaussian function as follows, such as 
that shown in Fig. 12(c) when the nonintegrability is intensified. 



P{Qm) = ^ , ^ (6-3) 



Where the parameter a is fitted by the general least squares procedure.^) 

In the case of the Henon-Heiles system at an energy of 0.1669, P{Qm) also 

shows a function in Eq. (6.3), as can be seen in Fig. 14(c). This corresponds to the 

fluctuation property Om in Fig. 13(c). 

Note that the statistical property of Om changes from a delta function to a 
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function in Eq. (6.3) with increasing perturbation strength. This highhghts the 
importance of Om as a fluctuating measure of non-integrabihty in the perturbed 
harmonic oscillator systems. 



Fig. 13. Computation of the fluctuation property Om = Rm x f m for the Henon-Heiles Hamiltonian 
with (a) E = 0.005, (b) E = 0.11796899, and (c) E = 0.1669. 



§7. Conclusions 

We conclude that this analysis provides us with a new method of detecting 
nonintegrabilities based on the degree of symmetry violation between and r^. 
It was demonstrated that slight changes in the properties of unstable orbits, which 
cannot be distinguished using the Lyapunov exponent, could be detected when using 
the absolute vahic of the correlation coefficient. 

We also verified numerically that the breaking of symmetry caused the prob- 
ability density for Qm to become the function in Eq. (6.3). It is hoped that a 
distribution function that describes the symmetry violation may be found in the 
future. Further research on the statistical principle that determines the probability 
distribution function for generating symmetry violation is needed; thus, statistical 
mechanics which can describe chaotic motion that breaks the symmetry based on 
this principle should be constructed. 
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Fig. 14. Probability density distribution of Qm- The width of the histogram dQm is 0.05. (a) 
E = -0.005, (b) E = 0.11796899, and (c) E = 0.1669. The dotted line is the function with the 
paremeter a = 9.9 in Eq. (6.3). 
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